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Summary 

Porcine epidemic diarrhoea virus (PEDV) and porcine deltacoronavirus (PDCoV) 
were first identified in Canada in 2014. Surveillance efforts have been instrumental 
in controlling both diseases. In this study, we provide an overview of surveillance 
components for the two diseases in Ontario (Canada), as well as PEDV and PDCoV 
incidence and prevalence measures. Swine herds located in the Province of Ontario, 
of any type, whose owners agreed to participate in a voluntary industry-led disease 
control programme (DCP) and with associated diagnostic or epidemiological informa¬ 
tion about the two swine coronaviruses, were eligible to be included for calculation 
of disease frequency at the provincial level. PEDV and PDCoV data stored in the 
industry DCP database were imported into the R statistical software and analysed 
to produce weekly frequency of incidence counts and prevalence counts, in addition 
to yearly herd-level incidence risk and prevalence between 2014 and 2016. The 
yearly herd-level incidence risk of PEDV, based on industry data, was 13.5%, 3.0% 
and 1.4% (95% Cl: 11.1-16.2, 2.0^4.2, 0.8-2.3), while the yearly herd-level inci¬ 
dence risk of PDCoV was 1.1%, 0.3%, and 0.1% (95% Cl: 0.5-2.2, 0.1-0.9, 0.0-0.5), 
for 2014, 2015 and 2016, respectively. Herd-level prevalence estimates for PEDV 
in the last week of 2014, 2015 and 2016 were 4.4%, 2.3% and 1.4%, respectively 
(95% Cl: 3.1-6.0, 1.5-3.3, 0.8-2.2), while herd-level prevalence estimates for PDCoV 
in the last week of 2014, 2015 and 2016 were 0.5%, 0.2% and 0.2%, respectively 
(95% Cl: 0.1-1.2, 0.0-0.6, 0.0-0.6). Collectively, our results point to low and 
decreasing incidence risk and prevalence for PEDV and PDCoV in Ontario, making 
both diseases possible candidates for disease elimination at the provincial level. 
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INTRODUCTION 


Porcine epidemic diarrhoea (PED) was first described in England in 
1971, and its causative agent, porcine epidemic diarrhoea virus 
(PEDV), was identified in 1978 (Chen et al., 2014). PEDV causes 
anorexia, vomiting, diarrhoea and dehydration in pigs, resulting in 


near 100% mortality for piglets during the first few days of life (Hill 
et al., 2014) and low mortality in older pigs. The virus spreads via 
the faecal-oral route, either through direct contact with an infected 
pig or through indirect contact with contaminated fomites. Wide¬ 
spread outbreaks were reported in Europe during the 1970s and 
1990s, while epidemics in Asia have caused significant disruption to 
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Asian pig production since 2008 (Williamson et al., 2013). Porcine 
deltacoronavirus (PDCoV, also known as swine deltacoronavirus, 
SDCV) was first identified in Hong Kong in 2012. The transmission 
modes and clinical signs due to PDCoV infections are similar to 
PEDV; however, the mortality rate is generally lower after PDCoV 
infections (Carvajal et al., 2015). 

PEDV emerged in North America in May 2013, while PDCoV 
was first confirmed in February 2014, both in the United States 
(Chen et al., 2014; Ma et al., 2015). These novel viruses rapidly 
disseminated throughout the US swine population, resulting in 
the mortality of an estimated seven million animals by May 2014 
(Jung & Saif, 2015). Due to the faecal-oral route of transmission, 
the infection spread through various mechanisms, including con¬ 
taminated transportation vehicles (Lowe et al., 2014). In Canada, 
PEDV emerged in January 2014 when a swine herd in Ontario 
tested positive for the virus (Kochhar, 2014). Imported spray- 
dried porcine plasma contaminated with PEDV was the likely 
pathway of introduction, as established through descriptive stud¬ 
ies (Pasma, Furness, Alves, & Aubry, 2016), analytical epidemio¬ 
logical studies (Aubry, Thompson, Pasma, Furness, & Tataryn, 
2017; O’Sullivan, 2015) and experimental investigations (Pasick 
et al., 2014). By July 2014, only 62 cases of PEDV had been 
reported in Ontario and the outbreak was largely under control 
(Pasma et al., 2016). 

Despite considerable impact on animal health in completely 
susceptible populations, the diseases caused by the two emerging 
porcine coronaviruses were not notifiable globally, according to 
the World Organization for Animal Health (OIE, 2017), and were 
not federally reportable or notifiable in Canada (Canadian Food 
Inspection Agency, 2014). At the provincial level, animal health 
regulations enabled consideration of PEDV and PDCoV as emerg¬ 
ing hazards (Government of Ontario, 2009), which allowed a mea¬ 
sured and appropriate response to the outbreak by the Ontario 
Ministry of Agriculture, Food and Rural Affairs (OMAFRA) via 
close collaboration with swine industry organizations. 

Initial emergence of the two novel porcine coronaviruses was 
followed by their successful elimination from several initial case 
farms (Misener, 2015). The high rate of successful PEDV elimination 
from individual herds, and effective measures that seemed to have 
minimized widespread viral dissemination, resulted in the current 
position of industry organizations that both infectious agents can 
and should be eliminated at the provincial level (Ontario Swine 
Health Advisory Board, 2017). 

A disease control programme (DCP) involves disease monitor¬ 
ing, surveillance, intervention and control strategies (Salman, 
2003). The DCP considered here has been voluntary in nature as 
defined elsewhere (Christensen, 2003). Furthermore, an important 
component of any disease control programme is measuring trends 
in incidence and prevalence, particularly when the disease of inter¬ 
est moves into the phase of possible elimination (Salman, 2003). 
With the infrastructure built for management of endemic diseases 
in Ontario, the data to support estimation of disease trends are 
available. Thus, the primary objective of this study was to 


estimate herd-level incidence and prevalence measures for PEDV 
and PDCoV in swine herds in Ontario (Canada) between January 
2014 and December 2016, based on industry data. The secondary 
objective was to describe relevant surveillance components that 
were used for identification of new PEDV cases. 
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MATERIALS AND METHODS 


2.1 


Data sources 


The Ontario Swine Health Advisory Board (OSHAB) maintains a 
database which contains premises information and PEDV/PDCoV 
herd status for producers enrolled in its voluntary regional disease 
control programmes. The program, and a database originally 
designed for management of porcine reproductive and respiratory 
syndrome virus (PRRSV) (Arruda, Poljak, Friendship, Carpenter, & 
Hand, 2015), was expanded to PEDV and PDCoV after their emer¬ 
gence. The data relevant for this work included unique identifiers, 
herd type, date of herd enrolment into the database, PEDV and 
PDCoV status of individual premises, and the date that individual 
premises changed their PEDV and PDCoV status. For inclusion into 
the study, swine herds could be of any herd type, but had to meet 
the following criteria: (i) be located in the Province of Ontario, (ii) 
participate in the voluntary industry-led disease control programme 
and (iii) have diagnostic or epidemiological information about the 
infection status of porcine epidemic diarrhoea virus (PEDV) or por¬ 
cine deltacoronavirus (PDCoV). The industry organization (i.e., 
OSHAB) provided relevant data to researchers for calculation of dis¬ 
ease frequencies under a separate data transfer agreement. 

Due to the voluntary nature of the DCP, the enrolment of herds 
into the database has been an ongoing process. This could have 
resulted in the date of enrolment being later than the date the infec¬ 
tion was originally detected in a specific herd. In rare instances 
where a herd’s enrolment date was not specified or occurred later 
than the first reported case of PEDV or PDCoV for that herd, the 
disease status date was entered as the herd enrolment date for the 
purposes of this report. This was done so that herd-level prevalence 
on a weekly basis could be properly calculated. An Open Database 
Connectivity (ODBC) connection to this database was established, 
and relevant tables were imported into R (R Core Team, 2016) using 
the rodbc library (Ripley & Lapsley, 2016). 


2.2 | Herd-level PEDV and PDCoV infection status 
definition 

The industry database maintained by OSHAB traces infection status 
of participating individual premises over time, which includes detec¬ 
tion of infection, as well as any subsequent infection, elimination or 
other change in infection status with the associated dates. In this 
database, premises can have four possible values for their infection 
status, as defined below: (i) confirmed positive —premises with herds 
that tested positive for PEDV/PDCoV based on laboratory test 
results from the Animal Health Laboratory (AHL) at the University of 
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Guelph (Guelph, ON, Canada). For a premise to be classified as con¬ 
firmed positive, it had to have an associated diagnostic submission 
number that includes at least one positive test for PEDV/PDCoV 
using RT-PCR, regardless of the number of specimens that were sub¬ 
mitted. An AHL reference number (also known as a G Number) was 
not available for some confirmed positive cases in the database, as 
the attending veterinarian obtained test results, but did not provide 
OSHAB with the AHL reference number; (ii) presumed positive —pre¬ 
mises which housed animals that were moved from positive sites at 
a prior stage in the production system (i.e., defined as positive due 
to pig movement). This information was obtained from attending 
veterinarians based on their knowledge of pig flow and movement, 
and was not based on diagnostic testing conducted on the premises 
of interest. For premises to be classified as presumed positive, the 
herd veterinarian simply needed to indicate that a specific site 
received pigs from PEDV/PDCoV-positive sites; (iii) presumed nega¬ 
tive —previously positive premises, either confirmed or presumed, 
that were tested using PCR tests according to industry guidelines 
and had all test results negative. The sampling requirements for 
declaring premises to be PEDV/PDCoV presumed negative varied 
based on the combination of herd type and the type of animal flow 
(i.e., all-in/all-out by barn, or continuous flow nursery and finisher 
herds). Complete criteria were, at the time of publishing, available on 
the website of the industry organization (Ontario Swine Health Advi¬ 
sory Board, 2015). Briefly, sampling strategy for farrow-wean, nurs¬ 
ery and finisher sites aimed to detect prevalence of virus-positive 
animals of at least 10%, with expected herd sensitivity of 95%, 
assumed test sensitivity of 98%, test specificity of 100%. In 
instances where the sampling material was oral fluid collected 
through cotton ropes, the assumption was that five pigs contributed 
oral fluids to one rope, and such fluids were considered a pooled 
sample. For farrow-finish or farrow-feeder sow sites, the same 
assumptions were made, except that the sampling strategy was 
required to detect prevalence of 5% with 95% herd sensitivity. 
Because of alternative strategies, the required sample size varied but 
a minimum sample size was four oral fluids. In addition, the testing 
time in sow herds was prescribed to be a minimum of 10 weeks 
post-infection and was required to be repeated three times in the 
case of farrow-wean (FW) sites, or two times in the case of farrow- 
finish (FF) or farrow-feeder (FG) sites. The recommended specimen 
type could be swab, Swiffer (for covering larger areas in a pen), or 
oral fluids, depending on the target age group. In FW herds, individ¬ 
ual farrowing crates were the target population for each individual 
sampling occasion, in particular if diarrhoea was evident. The mini¬ 
mum recommendation for one sampling occasion in FW herds was 
to sample four Swiffer samples, at least eight farrowing crates per 
one sample. Alternatively, individual swabs of 30 farrowing crates 
were deemed as acceptable sample after pooling 5:1. In FG farms, 
the recommended sample type was oral fluid, with recommendation 
to collect 12 oral fluids from nursery pigs. Similarly, in FF farms, six 
oral fluids were recommended for collection from nursery pigs and 
six for collection from finisher pigs. For all-in/all-out nursery and 


finisher farms, the recommendation was to sample six oral fluids, 
with added requirement that these herds should be supplied from 
sow herds with a confirmed negative status. All testing has been 
assumed to be performed using RT-PCR tests. Full description is 
available elsewhere (Ontario Swine Health Advisory Board (OSHAB), 
2015); (iv) confirmed negative —premises which have had no clinical 
signs or diagnostic evidence of PEDV/PDCoV for at least 6 months 
after the presumed negative status date. In addition, herds that were 
part of the Ontario voluntary DCP but were not tested for emerging 
porcine coronaviruses—due to lack of clinical or other types of diag¬ 
nostic or epidemiological triggers—had assigned status of NA (not 
available). 

2.3 | Weekly measures of disease frequency 

Weekly time series were then obtained for the three sets of statis¬ 
tics for both viruses: (i) number of premises with specific infection 
status in each week, (ii) number of new positive and new negative 
premises and (iii) number of premises in the database. 

For each week where a herd’s PEDV/PDCoV status was not 
reported, the status was set to the last-reported status using the zoo 
package (Zeileis & Grothendieck, 2005). For example, if a status is not 
reported for the current week, and “confirmed positive” was reported 
for the prior week, then the current week’s status is “confirmed posi¬ 
tive.” The individual premises data were then aggregated to counts of 
premises on a weekly basis. Based on the former time series, preva¬ 
lence area plots were generated, providing a visual assessment of 
“confirmed positive,” “presumed positive,” “presumed negative” and 
“confirmed negative” herds over time. The prevalence numerator was 
the positive herd count (sum of confirmed positive and presumed posi¬ 
tive herds) for a specific week, while the denominator was the herd 
count in the premises table for that week (calculated previously). 

Subsequently, the following disease status changes were tracked 
by week for each premise: (1) “Not Available” to “ Presumed Positive ,” 
(2) “Not Available" to “ Confirmed Positive,” (3) “Presumed Positive” to 
“Presumed Negative ,” (4) “Presumed Positive” to “Confirmed Negative,” 
(5) “Confirmed Positive” to “Presumed Negative,” (6) “Confirmed Posi¬ 
tive” to “Confirmed Negative,” (7) “Presumed Negative” to “Presumed 
Positive ,” (8) “Presumed Negative” to “Confirmed Positive,” (9) “Con¬ 
firmed Negative” to “Presumed Positive,” (10) “Confirmed Negative” to 
“Confirmed Positive.” 

Any status change leading to new presumed or confirmed positive 
status (i.e., status changes 1, 2, 7, 8, 9 and 10) was classified as new 
positive. Similarly, status changes 3, 4, 5 and 6 were classified as new 
negative. The number of new positives and new negatives was then 
aggregated to the weekly level throughout the study period. The latter 
time series were then used to construct a chart of cumulative inci¬ 
dence counts for each year, and epidemic curves were constructed for 
positive herds and herds which became negative. In addition, for each 
week, the incidence risk was calculated by dividing the number of 
cases that occurred in a specific week and by the number of herds that 
were eligible to become cases at the beginning of the week. 
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2.4 | Yearly measures of disease frequency 

Yearly herd-level incidence risk and prevalence at the end of the 
year were calculated in the following manner: the numerator for the 
yearly incidence risk calculation was defined as the cumulative num¬ 
ber of incident cases at the end of each year, while the denominator 
was defined as the number of herds in the DCP at the end of each 
year, minus the number of confirmed positive and presumed positive 
herds at the beginning of the specified year, minus half the number of 
herd additions to the database during the year (the number of addi¬ 
tions was treated as the number of withdrawals so this inherently 
open population could be converted to a closed population, and 
equations for calculation of risk used) (Eq. 1 (adapted from Dohoo, 
Martin, & Stryhn, 2003)): 

Yearly incidence risk (R) = N + PAR 


where N = cumulative incident cases at yearend and 


PAR (population at risk) = (number of herds in DCP at yearend) 

- (positive herds at beginning of year) 


— (- x number of additions to DCP 
during the year) 


( 1 ) 


For yearly herd-level prevalence, the numerator was defined as 
the total number of confirmed positive and presumed positive herds at 
the end of the year, while the denominator was defined as the num¬ 
ber of herds in the DCP at the end of each year. Both the yearly inci¬ 
dence risk and prevalence, originally calculated as proportions, were 
then multiplied by 100 to get the percentage values, with exact 95% 
confidence intervals obtained from the binom.test function in R. 

For the yearly incidence rate, we calculated total number of 
herd-years at risk for each year from the number of herds under risk 
in each week and used this as a denominator. The incidence rate 
was then expressed as number of cases per one herd-year. Exact 
95% confidence intervals on the incidence rate were obtained via 
the poisson.test function in R. 


2.5 | Surveillance components for PEDV 

Three surveillance components for identification of new cases were 
considered in this study. The cumulative number of presumed and 
confirmed cases between January 2014 and December 2016 main¬ 
tained in the OSHAB voluntary DCP database represents surveil¬ 
lance component 1. Surveillance component 2 is the official count of 
new cases maintained by OMAFRA and publicly available when this 
report was written (Ontario Pork, 2017). OMAFRA does not monitor 
PDCoV and defines a PEDV case as a herd which tests positive 
based on AHL’s RT-PCR test (Pasma et al., 2016), and is the first 
such case in the production system. In other words, the PEDV case 
count pertains to primary cases only; subsequent secondary cases 
due to animal movement in the production system are not included 
in OMAFRA reporting, although they could be confirmed as PEDV- 


positive in the diagnostic laboratory. Surveillance component 3 rep¬ 
resents cumulative number of positive submissions, due to any rea¬ 
son, maintained by the Animal Health Laboratory (AHL, University of 
Guelph), the largest diagnostic laboratory in Ontario for testing live¬ 
stock diseases. The overlap between the three surveillance compo¬ 
nents was calculated on the basis of demographic data. 
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RESULTS 


3.1 | Surveillance components for PEDV 


Figure 1 depicts three surveillance components considered in this 
study and their overlap. Briefly, during the 3-year period (January 
2014-December 2016), AHL reported a total of 974 positive sub¬ 
missions for PEDV (Fairies, personal communication 2017). Of the 
974 positive submissions, a cumulative total of 118 cases are 
reported in OMAFRA and OSHAB as new cases. The remaining 856 
positive submissions were not part of any surveillance component 
considered in this study. Of the 118 cases, 99 are captured by 
OMAFRA as primary PEDV cases; however, 70 of these 99 primary 
PEDV cases are also reported in the OSHAB database, leaving 29 
primary PEDV cases (99-70) reported in OMAFRA but not OSHAB. 



FIGURE 1 Venn diagram representing overlap among the three 
surveillance components for identification of new cases of the 
porcine epidemic diarrhoea virus (PEDV) in Ontario (Canada) 
between 2014 and 2016—The Animal Health Laboratory (AHL), 
Ontario Ministry of Agriculture Food and Rural Affairs (OMAFRA) 
and the Ontario Swine Health Advisory Board (OSHAB). The size of 
individual circles should be proportional to the number of positive 
cases recorded by a specific surveillance component. The 
overlapping areas should represent cases that were identified by as 
many surveillance components as are being overlapped. The 
numbers listed in a specific area represent number of cases that are 
identified by a specific surveillance component or combination 
thereof 
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TABLE 1 Herd-level incidence risk and rate of two novel porcine coronaviruses (PEDV and PDCoV) in Ontario swine herds between 2014 
and 2016, and estimated prevalence of positive cases at the end of each year based on data provided in the Ontario Swine Health Advisory 
Board (OSHAB) Disease Control Program (DCP) database (average number of herds for 2014-2016 = 1093) 


Year 

Cumulative n 

of new cases 

Incidence 
risk (%) 

95% Cl (%) 

Incidence rate (cases 95% Cl (cases 
per herd-year) per herd-year) 

Number of cases Prevalence at 
at year-end year-end (%) 

95% Cl 

(%) 

Porcine epidemic diarrhoea virus 

2014 

95 

13.49 

(11.06-16.24) 

0.14 

(0.12-0.18) 

36 

4.36 

(3.07-5.99) 

2015 

29 

2.97 

(2.00^1.24) 

0.03 

(0.02-0.05) 

27 

2.25 

(1.49-3.26) 

2016 

17 

1.42 

(0.83-2.26) 

0.02 

(0.01-0.03) 

17 

1.35 

(0.79-2.16) 

Porcine deltacoronavirus 

2014 

8 

1.14 

(0.49-2.23) 

0.011 

(0.005-0.022) 

4 

0.48 

(0.13-1.24) 

2015 

3 

0.30 

(0.06-0.87) 

0.003 

(0.001-0.009) 

2 

0.17 

(0.02-0.60) 

2016 

1 

0.08 

(0.00-0.45) 

0.001 

(0.000-0.005) 

2 

0.16 

(0.02-0.57) 


The remaining 19 cases (974-856-99) are non-primary PEDV-posi- 
tive cases (i.e., not the first outbreak in the associated production 
system); hence, they are captured in AHL and reported in the 
OSHAB voluntary DCP database, but not OMAFRA. Lastly, 52 
PEDV-positive cases are strictly in OSHAB and are not reported in 
any other surveillance component. It should be noted that these 52 
cases are “presumed positive” due to pig movement, and as such 
have no supporting AHL test result, an attribute which also makes 
them ineligible for OMAFRA reporting. The OSHAB total for PEDV 
cases is 52 + 19 + 70 = 141, which is also the cumulative number 
of incident cases over the 3-year period. 


could be observed that PEDV had a rapid increase in the number of 
prevalent cases until March 2014, which then peaked in the summer 
of the same year, only to be followed by a decline. Figure 5 displays 
the weekly development of infection status as a proportion. Note 
that in the latter figure, the denominator changed on a weekly basis 
as the number of herds in the voluntary DCP changed. Figures 6 
and 7 display the number of premises in a distinct PDCoV status as 
a count (Figure 6) and proportion (Figure 7). The weekly incidence 
risk of PEDV and PDCoV in 2014—2016 is provided in the supple¬ 
mentary material (Figures SI and S2, respectively). 


3.2 | Yearly incidence and prevalence measures 


4 


DISCUSSION 


Table 1 contains number of new infected premises for the two 
viruses in each year based on the industry’s voluntary DCP data¬ 
base, number of positive cases in the last week of each year, as 
well as incidence risk and prevalence estimates. In addition, 95% 
confidence intervals are provided, recognizing that this surveillance 
component was a census of herds participating in the voluntary 
DCP, but does not contain all herds from the source population. In 
brief, annual incidence risk was the highest in 2014 for PEDV 
(13.49%) and PDCoV (1.14%), and the lowest in 2016 for PEDV 
(1.42%) and PDCoV (0.08%). The end of year prevalence for PEDV 
ranged between 4.4% and 1.4% in 2014 and 2016, respectively. A 
similar trend can be seen for PDCoV (Table 1). Furthermore, similar 
values have been observed for incidence rates in each year 
(Table 1). 

3.3 | Weekly incidence and prevalence measures 

The cumulative number of new PEDV and PDCoV cases per week 
for the 3 years is provided in Figures 2 and 3, respectively. Notably, 
there were many new PEDV cases detected in the winter of 2014 
when the disease was first introduced to Ontario and Canada. 

Figure 4 depicts the number of herds categorized into one of the 
four groups with respect to PEDV infection. From this figure, it 


Emerging infections considered in this work are examples of produc¬ 
tion-limiting diseases that are not federally reportable (Canadian 
Food Inspection Agency, 2017) or listed in an OIE list of diseases 
(World Organization for Animal Health, 2017). Consequently, 
resources available through response by the regulatory veterinary 
authorities could be limited. In the case of the Ontario outbreak, the 
regulatory framework existed at the provincial level because PEDV 
was considered an emerging hazard under the provincial Animal 
Health Act (Government of Ontario, 2009) and because of suspicion 
of feed involvement, which is under the jurisdiction of national vet¬ 
erinary authorities. The response to the outbreak was collaborative 
in nature and is best depicted by the evaluation of surveillance com¬ 
ponents for PEDV. 

The OMAFRA and industry approaches were different with 
respect to surveillance coverage, case definition and follow-up 
time. OMAFRA surveillance was mostly concerned about identifi¬ 
cation of new cases and their investigation in the early phase of 
the outbreak. It contained a census of all PEDV cases, and all 
herds from the source population were eligible to be listed as a 
case, at least during the initial phase of the outbreak when PEDV 
could be considered an emerging hazard. However, only the first 
case in a given production system was counted as a case and 
secondary cases due to planned animal movement to other 












AJAYI ET AL. 


6 


Wiley- 


Iransboundary and Emerging Diseases 



FIGURE 2 Cumulative weekly confirmed and presumed positive cases for porcine epidemic diarrhoea virus (PEDV) in Ontario (Canada) for 
2014-2016, based on data provided in the Ontario Swine Health Advisory Board (OSHAB) Disease Control Program (DCP) database (average 
number of herds for 2014—2016 = 1093) 



FIGURE 3 Cumulative weekly confirmed and presumed positive cases for porcine deltacoronavirus (PDCoV) in Ontario (Canada) for 2014— 
2016, based on data provided in the Ontario Swine Health Advisory Board (OSHAB) Disease Control Program (DCP) database (average 
number of herds for 2014-2016 = 1093) 


premises were not counted as additional cases. For example, if 
PEDV was detected in a sow herd in a multisite production 
system, the sow herd would be counted as a case, but not the 
nursery or finisher sites supplied from the sow herd. For this 
surveillance system, all cases had to have diagnostic confirmation 


through laboratory submissions. For this reason, the OMAFRA 
PEDV surveillance system is completely nested within the diagnos¬ 
tic laboratory. It is notable that -17.1% of PEDV cases were iden¬ 
tified through this surveillance system only, which is reflective of 
the reportable nature (provincially) of the emerging hazard. 
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FIGURE 4 Number of Herds with confirmed positive, presumed positive, presumed negative or confirmed negative status for porcine 
epidemic diarrhoea virus (PEDV) in Ontario (Canada) for 2014-2016, based on data provided in the Ontario Swine Health Advisory Board 
(OSHAB) Disease Control Program (DCP) database (average number of herds for 2014-2016 = 1093) 
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FIGURE 5 Percentage of Herds with confirmed positive, presumed positive, presumed negative or confirmed negative status for porcine 
epidemic diarrhoea virus (PEDV) in Ontario (Canada) for 2014-2016, based on data provided in the Ontario Swine Health Advisory Board 
(OSHAB) Disease Control Program (DCP) database (average number of herds for 2014—2016 = 1093) 


The OSHAB surveillance system was concerned with identifica¬ 
tion of new cases, reducing dissemination of infection through ani¬ 
mal movement, and following the infection status over time for the 
purposes of monitoring disease trends over time. As a voluntary pro¬ 
gramme, the source population consisted of premises that volun¬ 
teered to participate. At the end of 2016, the number of sites in the 
OSHAB database was 1255, which represents approximately 49% of 


2,556 active swine sites in Ontario (Brisson, 2014). Thus, the 
OSHAB surveillance coverage was lower than for the OMAFRA 
surveillance, which had 100% coverage—by law, PEDV-infected 
herds were reported to OMAFRA during the phase when the hazard 
was still considered as emerging. 

The OSHAB case definition also included secondary sites that 
were confirmed positive due to animal movement and had 

































AJAYI ET AL. 


8 


Wiley- 


fransboundary and Emerging Diseases 


18- 

16- 



2014 2014 2015 2015 2016 2016 


Confirmed Positive Presumed Positive Presumed Negative Confirmed Negative 

FIGURE 6 Number of Herds with confirmed positive, presumed positive, presumed negative or confirmed negative status for porcine 
deltacoronavirus (PDCoV) in Ontario (Canada) for 2014—2016, based on data provided in the Ontario Swine Health Advisory Board (OSHAB) 
Disease Control Program (DCP) database (average number of herds for 2014-2016 = 1093) 
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FIGURE 7 Percentage of Herds with confirmed positive, presumed positive, presumed negative or confirmed negative status for porcine 
deltacoronavirus (PDCoV) in Ontario (Canada) for 2014—2016, based on data provided in the Ontario Swine Health Advisory Board (OSHAB) 
Disease Control Program (DCP) database (average number of herds for 2014-2016 = 1093) 


laboratory confirmation, as well as secondary sites where movement 
of PEDV-positive animals occurred although testing was not carried 
out This approach allowed identification of -41.8% of cases that 
would otherwise be unaccounted for in the OMAFRA system 
(n = 71), which included 30.6% of total cases (n = 52) that were not 
confirmed through diagnostic testing (Figure 1). Given the purpose 
of the industry’s surveillance system, recording infection in the 


secondary sites was a logical choice as they could still contribute to 
disease dissemination through multiple pathways. Such premises had 
to be declared as presumed or confirmed negative, based on diag¬ 
nostic testing that met pre-determined criteria conducted on the 
premises. 

The Venn diagram of surveillance components (Figure 1) also 
shows a large volume of positive test results for PEDV in a 
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laboratory that did not contribute to the identification of new cases. 
Such testing serves different purposes and likely contains additional 
surveillance components which are not considered here. An example 
is testing trailers in abattoirs for the purposes of case identification 
through risk-based approaches. Surveillance systems for emerging 
diseases could benefit from a clearer definition of such surveillance 
components, so that diagnostic tests could be easily aggregated into 
appropriate streams. 

The collaborative efforts at controlling PEDV are most clearly 
seen in the incidence counts for 2014. As expected at the onset of 
any outbreak, particularly due to massive contamination through a 
common source, there was a rapid climb in the number of infected 
farms. However, this outbreak was brought under control relatively 
quickly. By May 2014, the worst of the outbreak was over, and 
there were only single digit weekly infections (no more than 5) for 
the remainder of the year. In addition, the initially positive herds 
which contributed to the rapid climb in infections between January 
and May 2014 became negative within a short period and con¬ 
tributed to the rapid increase in negative herds between June and 
December 2014. 

The incidence counts for 2015 reveal a more proactive approach 
to OSHAB PEDV surveillance and disease control. Not only were 
there substantially fewer infections (infections were only reported 
for 15 weeks out of the year, and there were at most 5 infections 
per week), the number of negative herds also kept up with the num¬ 
ber of positive herds (for 9 of the 15 weeks where infections were 
reported, there were also negative herds). Overall, the number of 
negative herds in 2015 outpaced the number of positive herds, indi¬ 
cating proactive steps to minimize the number of new infections 
while aggressively taking steps to eradicate existing ones. This 
stands in stark contrast to 2014, where disease eradication efforts 
only seemed to take on momentum towards the latter half of the 
year. 

The trend in incidence counts from 2015 continued into 2016, 
with lower numbers for positive herds and the number of negative 
herds outpacing positive ones. However, in contrast to 2015, 
roughly 75% of new infections occurred in the first quarter of the 
year (January-March), while 80% of new negative herds were 
reported in the last quarter of the year (late September-December). 
The cluster of negative herds at the end of the year suggests that 
with new infections under control, emphasis was placed on the man¬ 
agement of prevalent cases, taking into consideration lower viral 
transmission during warmer periods of the year. 

The prevalence plots reveal ebbs and inclines in PEDV preva¬ 
lence, in what appears to be a cyclical pattern: a peak in July 2014, 
a low in January 2015, a peak in July 2015, a low in October 2015 
and a peak in February 2016 (Figure 4). The numbers were relatively 
constant from late February to July 2016, declined in late July and 
slowly declined until the end of the year. Given the prevalence fluc¬ 
tuations between 1% and 2% for September 2015-December 2016 
(Figure 5), and the repeated patterns, there is a danger that PED 
may become an endemic disease with low prevalence and limited 
pockets of infection in Ontario. As such, further investigation and 


action are needed to ensure that the peaks are replaced by a steady 
decline. 

As the same disease prevention and control measures for PEDV 
also apply to PDCoV (Ontario Ministry of Agriculture Food and Rural 
Affairs, 2014), there were few reported infections for PDCoV (essen¬ 
tially, controlling for PEDV also controls for PDCoV). It is therefore 
not surprising that the incidence count patterns are similar: in 2014, 
there were more PDCoV-positive herds than negative herds, while in 
2015, there were more negative herds than positive ones. However, 
in 2016 there was just one new PDCoV infection, with a single herd 
also becoming negative, an indication that if current disease control 
efforts continue, there may be no new infections in subsequent 
years. 

The prevalence figures for PDCoV are quite different from 
PEDV—rather than a cyclical pattern as previously observed, the 
numbers are more erratic, with sudden peaks followed by con¬ 
stants over prolonged periods, and then declines. The relatively 
static and low prevalence (0.1%-0.2%) observed between July 
2015 and December 2016 suggests that PDCoV is a candidate 
for disease elimination. There is also a possibility that PDCoV 
cases are underreported by the industry, either because of poten¬ 
tially lower clinical impact in swine herds or because of perceived 
lower importance than PEDV. Nonetheless, when any of the three 
porcine coronaviruses (PEDV, PDCoV, TGEV—transmissible gas¬ 
troenteritis virus) is suspected in a herd and diagnostic material is 
submitted to the AHL, the diagnostic testing is automatically con¬ 
ducted for all three viruses. 

Some limitations of the current study include the fact that 
PEDV and PDCoV cases are actively pursued for inclusion into the 
voluntary DCP by industry organizations. Such strategy is likely to 
result in estimates of incidence and prevalence measures that are 
higher than in the source population. Because of the inherently 
open nature of the voluntary DCP, we had to modify the formulae 
for calculation of incidence risk. It could be argued that presumed 
cases are not diagnostically confirmed and are therefore subject to 
misclassification. However, the reality of a voluntary DCP for pro¬ 
duction-limiting diseases is that resources to conduct large-scale 
testing are scarce and, as such, need to be carefully deployed. Fur¬ 
thermore, in order to confirm premises as presumed negative, diag¬ 
nostic testing to confirm absence of infection at the design 
prevalence level is still required. Also, the criteria to declare con¬ 
firmed negative status is arguably open-ended and could be further 
improved. 

In conclusion, this study provides estimates of incidence and 
prevalence measures in Ontario based on industry data collected 
through voluntary disease control programmes. The data suggest 
that annual incidence risk and prevalence estimates are low and 
have been steadily decreasing between 2014 and 2016 for PEDV 
and PDCoV. Current estimates of disease frequency support plan¬ 
ning of disease elimination at the provincial level, but much informa¬ 
tion should be available about factors that led to time to elimination 
in individual herds. In addition, our evaluation of surveillance compo¬ 
nents indicates that the two surveillance components were 
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complementary and focused on different aspects of surveillance. 
OMAFRA surveillance was mostly focused on identification of pri¬ 
mary cases aimed at quick disease investigations and traceability in 
the face of the outbreak, whereas OSHAB surveillance has the 
added benefit of having sufficient data that allow long-term evalua¬ 
tion of disease trends, long-term disease management and tracing 
disease status of individual herds over time. The OSHAB voluntary 
DCP database also provides a good tool for calculating weekly 
prevalence and incidence measures, which is a valuable statistic for 
producers and animal health experts during all phases of disease out¬ 
break and control. 
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